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ABSTRACT 



We have discovered 12 new molecular hydrogen emission-line objects (MHOs) 
in the vicinity of the candidate massive young stellar object Mol 121, in addition 
to five that were previously known. H2 2.12-/xm/H2 2.25-//m fiux ratios indi- 
cate another region dominated by fluorescence from a photo-dissociation region 
(PDR), and one region that displays an anomalously low H2 2.12-/im/H2 2.25-/im 
flux ratio (<1) and coincides with a previously reported deeply embedded source 
(DES). Continuum observations at 3 mm reveal flve dense cores; the brightest 
core is coincident with the DES. The next brightest cores are both associated 
with cm continuum emission. One of these is coincident with the IRAS source; 
the other lies at the centroid of a compact outflow defined by bipolar MHOs. 
The brighter of these bipolar MHOs exhibits [Fe II] emission and both MHOs 
are associated with CH3OH maser emission observed at 95 GHz and 44 GHz. 
Masses and column densities of all five cores are consistent with theoretical pre- 
dictions for massive star formation. Although it is impossible to associate all 
MHOs with driving sources in this region, it is evident that there are several 
outflows along difi'erent position angles, and some unambiguous associations can 
be made. We discuss implications of observed H2 2.12-//m/H2 2.25-//m and [Fe II] 
1.64-//m/H2 2.12-//m flux ratios and compare the estimated total H2 luminosity 
with the bolometric luminosity of the region. We conclude that the outflows are 
driven by massive young stellar objects embedded in cores that are likely to be 
in different evolutionary stages. 

Subject headings: ISM: individual objects (Mol 121/ IRAS 20188-^3928/ 
G077.4622-I-01.7600) — ISM: jets and outflows — stars: formation — stars: massive — 
stars: pre-main sequence — stars: protostars 
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1. INTRODUCTION 

Outflows play a critical role in the formation of isolated low-mass stars throughout the 
embedded accretion phases of evolution. They are thought to be driven by magnetic stresses 
in accretion disks that eject some of the inflowing matter and carry off angular momentum 
(e.g., Shu et al. 1994; Pudritz et al. 2006; Matsakos et al. 2009). In contrast, the role of 
outflows in the formation of stars > 8 Mq remains elusive, because these stars are rare, 
form in tight clusters along with low- and intermediate-mass stars that may be in different 
stages of evolution, lie in heavily obscured regions often at distances greater than several 
kiloparsecs, and may be evolving on dynamical time scales less than a few thousand years 
(e.g., Zinnecker & Yorke 2007). Large surveys with single-antenna microwave telescopes 
indicate that energetic molecular outflows are present in massive star-forming regions (e.g., 
Zhang et al. 2005), but it is not generally possible to link individual outflows with driving 
sources, even with the current generation of millimeter- wave interferometers. 

The situation is somewhat better in the near-infrared, where the H2 1-0 S(l) line at 
2.12 iim is a particularly powerful tracer of shocks in molecular outflows (Davis et al. 2010 
and references therein; Smith 2012). The current generation of sensitive wide-field near 
infrared imagers enables large surveys of massive star-forming regions to be conducted at 
arcsecond or sub-arcsecond resolution in reasonable periods of observing time. Using the 
Apache Point Observatory (APO) 3.5-m telescope, Wolf-Chase et al. (2012, in preparation) 
have conducted a narrow-band near-infrared survey for shocked emission in 28 massive 
star-forming regions that exhibit energetic molecular outflows (Zhang et al. 2005). In 
Wolf-Chase et al. (2012, hereafter WSS12), we presented near-infrared results for Mol 160 
(IRAS 23385-1-6053), along with 3-mm continuum and line observations obtained with the 
Combined Array for Research in Millimeter- wave Astronomy (CARMA). In this paper, 
we present narrowband, ~ l"-resolution near-infrared APO images and ~ 2"-resolution 
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CARMA 3-mm continuum and line observations for Mol 121 20188+3928). 

Mol 121 has a spectral energy distribution (SED) resembling an ultra-compact (UC) 
H II region and is associated with radio continuum emission (Wood & Churchwell 1989; 
Molinari et al. 1998). Varricatt et al. (2010) included it in a NIR imaging survey of 50 
regions containing intermediate- and high-mass YSO outflow candidates. They present an 
overview of background observations of this object in the appendix to their paper (A37) and 
note that prior distance estimates to Mol 121 span an order of magnitude (0.3 - 3.91 kpc). 
Mol 121 (MSX6C designation G077.4622+01.7600) is categorized in the Red MSX Survey 
(RMS) as Class H II/YSO, and is among the RMS objects that have uniquely constrained 
kinematic distances and good SED fits (Mottram et al. 2011, Table 1). In this paper, we 
assume the RMS catalog distance of 1.7 kpc and bolometric luminosity of 8.2x10^ Lq. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Apache Point Observatory: NICFPS 

We obtained narrowband, near-infrared images of Mol 121 using the Near-Infrared 
Camera and Fabry-Perot Spectrometer (NICFPS) on the Astrophysical Research 
Consortium (ARC) 3.5-m telescope at the APO in Sunspot, NM on UT 2008 Oct 13. We 
used the following filters and total integration times: H2 2.12-yum (1740 s), H2 2.25-/im 
(1740 s), and [Fe II] 1.64-/xm (1800 s). We also used two narrowband filters to allow for 
continuum subtraction in the final images: an H2-continuum filter centered on 2.13 /im 
(1740 s) and an [Fe Il]-continuum filter centered on 1.65 //m (1800 s). Our basic data 
acquisition, reduction, and calibration methods are described in WSS12. 

Thirteen isolated stars having the best 2MASS photometry fiags were chosen in the 
narrow-band line images as well as in the corresponding narrow-band continuum images. 
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We performed aperture photometry on these stars using the DAOPHOT package in IRAF 
with an aperture radius of 12 pixels, ~4 times the FWHM, since the seeing was ~1" in each 
image. The continuum images were scaled using a scaling factor and then subtracted from 
the narrow-band line images so that on average the stars disappear. The scaling factor was 
chosen so that the fractional residual fluxes of the stars within the aperture are as close 
to zero as possible in the resulting continuum-subtracted image. This procedure worked 
well except for the brightest stars in the field, which show small residual artifacts. We flux 
calibrated our data using equations (1) & (2) from WSS12 with dX = 5.5 nm for the [Fe II] 
1.64-/im filter (Apache Point Observatory 2010). 

Since the H2 sources can have irregular morphologies, we performed irregular aperture 
photometry using the highly flexible kan^ software version 1.3. The apertures were 
defined using polygons. For each aperture in each image, we drew ten background regions. 
These background regions were chosen so they sample the background as close to the 
aperture as possible, and were roughly of the same size as the apertures. We then used 
the average background-subtracted flux within each aperture for subsequent calculations. 
The uncertainties in Tables 1 & 2 correspond to ±lcr, and are the combinations of 
three categories of uncertainty; calibration, subtraction, and aperture. The calibration 
uncertainty incorporates the 2MASS catalog K,, magnitude uncertainties, our instrumental 
stellar magnitude uncertainties and the standard deviation for the zero-point correction. 
The subtraction uncertainty is estimated using the standard deviation from zero for the 
fractional residual fluxes in the continuum subtraction procedure. The aperture uncertainty 
is the combination of photon noise within the aperture and background regions and the 
standard deviation from the ten background regions used for each aperture. 



"^http:/ /www. bu.edu/iar/kang/ 
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2.2. CARMA 

We observed a field around Mol 121 approximately 1.5' across with CARMA on UT 
2010 March 11. The interferometer was in the C configuration with baselines ranging from 
26-370 meters. We simultaneously observed the continuum and the CH3OH J=8^7 line 
at 95.169 GHz. Two spectral bands were configured for maximum continuum bandwidth, 
covering 468 MHz each, for a total of 1.872 GHz (including both sidebands), on either side 
of a central frequency of 93.585 GHz. The band observing CH3OH was in the upper side 
band with a 31-MHz bandwidth and 63 channels, resulting in a velocity resolution of 1.538 
km s~^. The passband calibrator was 1927-1-730, the fiux calibrator was MWG349, and the 
phase calibrator was 2007-1-404. The total integration time on source was 4.92 hours. The 
size (FWHM) of the synthesized beam was 1.967" x 1.516". We reduced and calibrated the 
data as described in WSS12. 

3. RESULTS 

3.1. H2 & [Fe II] Emission 

Figures 1-3 present our narrow-band NIR images. Our H2 2.12-//m, H2 2.25-//m, 
[Fe II], and continuum-subtracted line images are shown in Figure 1. Figure 2 displays 
the continuum-subtracted H2 2.12-/xm image with the regions used to identify molecular 
hydrogen emission-hne objects (MHOs), compact emission-hne sources thought to be 
associated with outfiows. Two other emission features that are not MHOs are also 
identified: fluorescent gas in a photo-dissociation region (PDR) and a region where the 
H2 2.12 /im/2.25 /im fiux ratio is less than one (marked as 'Anomaly'). The anomalous 
H2 emission is discussed further in §4.1. For the MHOs, we follow the numbering scheme 
of the on-hne catalog of these objects, hosted by the Joint Astronomy Gentre in Hawaii 
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(Davis et al. 2010). In all, we identify 12 new MHOs in addition to five that have been 
catalogued previously (see Figure 1: bottom left panel, Figure 2, & Table 1). The feature 
originally identified as MHO 867 by Varricatt et al. (2010) actually encompasses both the 
PDR and the bow-shaped feature we associate with the MHO. Although the IRAS source 
appears to be offset at the eastern edge of the PDR in the H2 2.12-//m and [Fe II] images, it 
is centrally located in the H2 2.25-/im images (Figures 1: top and bottom middle panels) 
and the 3-color Wide-field Infrared Survey Explorer (WISE: Wright et al. 2010) image 
presented in Figure 4. The presence of shocked gas to the west of the IRAS source may 
explain the enhancement of H2 2.12-//m and [Fe II] emission in this direction. 

Table 1 lists the MHO designations (column 1), areas used in the aperture photometry 
(column 2), positions of the peak emission (columns 3 & 4), H2 2.12-/im line fluxes (column 
5), H2 2.25-/xm line fluxes (column 6), and H2 2.12-/im/H2 2.25-/xm flux ratios (column 7). 
Columns 8—10 present results for fluxes integrated over 2-pixel circular apertures centered 
on 2.25-/im peaks. Where the H2 2.25-//m line flux is below 3 times the noise {3a) within 
an aperture, the H2 2.25-/xm line is regarded as undetected and no value is given in columns 
6 or 9. In these cases, lower limits to the H2 2.12-/xm/H2 2.25-/xm flux ratios are calculated 
using the 3a value for that aperture as an upper limit for the H2 2.25-//m line flux. The 
difference in extinction at 2.12 /xm and 2.25 /xm is very small, only a tenth of a magnitude 
(assuming ~ A~^'^), yielding a corresponding brightness ratio of ~1.1, which has a 
neghgible effect on the flux ratios. 

The expected H2 2.12-/im/2.25-/xm line ratio for UV excitation (as in PDRs) is ~ 1.9. 
It is significantly higher for H2 excited by shocks, with a ratio of ~ 3 - 20 for different types 
of shocks (e.g.. Black & van Dishoeck 1987; Smith 1994; Gredel & Dalgarno 1995; Smith 
1995; Smith et al. 2003). Two types of shocks are observed in outflows from star-forming 
regions. Continuous (or C-type) shocks occur when the shock is magnetically cushioned 
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by a relatively high transverse magnetic field, such that the shock thickness is broad and 
deceleration of gas is relatively gradual, resulting in lower excitation and temperatures than 
in Jump (or J-type) shocks. In the latter, the gas undergoes a sharp increase in temperature 
within a narrow zone, where molecules are dissociated (e.g., Barsony et al. 2010; Smith 
2012). 

The observed line ratios for objects listed as MHOs are consistent with shock excitation, 
but not for the PDR and Anomaly. Since the high flux ratios for many of the MHOs may 
reflect the fact that 2.25-/xm emission is not detected over the full extent of the apertures 
used for the 2.12-/xm emission (see Table 1, columns 5-7), we also computed flux ratios 
using a circular aperture of radius 2 pixels centered on the 2.25-//m emission peaks (see 
Table 1, columns 8-10). A comparison of columns 10 and 7 indicates that although H2 
2.12-/xm/2.25-/xm flux ratios are typically lower in the 2-pixel radius apertures, the ratios 
are still consistent with shock excitation. 

The [Fe II] 1.64-/xm line probes fast, dissociative shocks with velocities > 50 km s~^. 
When it is observed in jets from low- and intermediate-luminosity YSOs, it is typically 
localized in compact spots, representing only a small fraction of the region emitting in H2, 
and its presence docs not appear to correlate strictly with the evolutionary stage of the 
driving source (Caratti o Garatti et al. 2006). Fewer [Fe II] observations have been reported 
for high-luminosity objects, in part due to the greater extinction at 1.64 /im. Curiously, 
Caratti o Garatti et al. (2008) found that the only [Fe II] emission toward the H2 jet 
associated with the Massive Young Stellar Object (MYSO) IRAS 20126+4104 is close to 
the source. Table 2 lists the three [Fe II] emission-hne features we detected. The flrst two 
features are associated with MHO 864, and the third feature is associated with, and largely 
overlaps, MHO 867, the bright bow- or 'bean'-shapcd feature that opens toward the IRAS 
source. In MHO 864, the two [Fe II] emission-line features lie along the inferred jet axis of 
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the outflow. One is located at the apex of the MHO and the other is close to core D, which 
presumably harbors the driving source of the outflow (see Figure 2 and §4.2.4). In MHO 
867, the morphology of the [Fe II] emission is similar to the MHO, but over a smaller region. 
To investigate excitation conditions in these MHOs, we computed line ratios in 2-pixel 
radius apertures along the major and minor axes of MHO 864, and along and transverse to 
the bow-shaped structure of MHO 867. Placement of these apertures is indicated in Figure 
3 and emission line flux ratios are presented in Table 3. For the apertures with H2 2.25-/xm 
detections, 2.12 //m/2.25 iim flux ratios range from 6.8 - 14.5. 1.64-//m/2.12-//m flux ratios 
range from 0.03 - 0.9 for MHO 864 apertures and 0.5 - 1.9 for MHO 867 apertures. We 
interpret these results in §4.1. 

The terms "Green Fuzzies" and "Extended Green Objects" (EGOs) have been apphed 
to regions of extended 4.5-/xm emission that are often associated with shocked H2 emission 
from outflows (e.g., Cyganowski et al. 2008). The nomenclature reflects the use of green 
as the representative color for the InfraRed Array Camera (IRAC: Fazio et al. 2004) 
4.5-/im band in Spitzer Space Telescope 3-color images produced from Galactic Legacy 
Infrared Mid-Plane Survey Extraordinaire observations (GLIMPSE: Benjamin et al. 2003; 
Churchwell et al. 2009). Figure 4 uses a similar 3-color scheme to represent the WISE 3.4, 
4.6, & 12-/xm bands. Here, objects that are bright in the 4.6-/xm band appear to be green. 
MHOs 864, 865, & 866 are the only MHOs that are clearly prominent in the 4.6-/xm band. 
MHOs 864 & 865 have the highest H2 2.12-//m fluxes in Mol 121. On the other hand, there 
are several patches of faint 4.6-//m emission that do not correspond to MHOs. We note that 
this result is similar to that of Lee et al. (2012), who imaged 34 Spitzer GLIMPSE EGOs 
with the United Kingdom Infrared Telescope (UKIRT) and noted that shocked molecular 
hydrogen emission is typically more extended than EGOs and that EGOs occasionally 
trace scattered continuum hght (see also Barsony et al. 2010). This suggests that although 
the 4.5-/xm & 4.6-/xm bands may be useful in identifying features associated with bright 



outflows, they are not always dominated by shocked outflow emission. 



3.2. 3-mm Continuum Emission 

Our CARMA 3-mm continuum data reveal five cores in the Mol 121 region. Figure 5 
shows contours for the 3-mm continuum emission overlaid on the relevant portion of our 
continuum-subtracted H2 2.12-//m image. Table 4 lists physical properties of the cores. 
Assuming a distance of 1.7 kpc, deconvolved core diameters, {DjnaxDmin)K range from 790 
AU - 3500 AU. The cores are discussed individually in §4.2. 



3.2.1. Estimating Core Mass & Column Density 

Assuming the dust emission is optically thin at 3 mm, the mass of a core is calculated 
in the following manner: 



M ^'^^^ 



Here, d is the distance to the source, Sj, is the total flux density, R is the gas-to-dust 
ratio, B^{Td) is the Planck function at dust temperature Tp, and Ki, is the dust opacity 
per gram of gas (e.g., Enoch et al. 2007). We use k.^ = 0.178 cm^ g~^, the value obtained 
by extrapolating from Table 1 in Ossenkopf & Henning (1994) using a power law {k oc A~^) 
with (3 — 1.8, and assuming dust grains with thin ice mantles coagulated for 10^ yr at a 
gas density of 10^ cm~^. Typically, dust temperatures in protostellar cores are assumed 
to be ~ 20 K (e.g., Rosolowsky et al. 2010; Lee et al. 2011; Chen et al. 2012), but dust 
temperatures in PDRs and some massive protostellar cores are closer to 30 K (e.g., Molinari 
et al. 2008; Anderson et al. 2012), and Enoch et al. (2007) note that plausible values might 
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vary from 5-30 K. The SED of IRAS 20188+3928 has been fit with a cold dust component 
of 39 K and warm dust component of 100 K (McCutcheon et al. 1995). Most of the cores 
detected by CARMA are probably colder than this, so we have opted to calculated masses 
for = 20 K & 40 K (Table 3, columns 7 & 8). For T^ = 40 K, derived masses are ~ 2.1 
times lower than for Td — 20 K. 

We calculate the column density from the peak flux density Ii, in a similar manner to 
Enoch et al. (2007): 



N{H2) = ^ (2) 

Here, R, k^, and B^,{T£))-, are the same as in equation (1), Vtmh is the size of CARMA's 
main beam at 3 mm, and m^/j is the mass of molecular hydrogen. Columns 9 & 10 in 
Table 3 show the results of this calculation for both dust temperatures used in the mass 
calculation. 

Theoretical work predicts a mass column density threshold > 1 g cm~^ for massive star 
formation (Krumholz & McKee 2008; Krumholz et al. 2010). The mass column density is 
the column density N{H2) multiplied by 111^2 and a factor of 1.36 to account for elements 
heavier than hydrogen (Simon et al. 2001). Regions where intermediate mass stars are 
thought to be forming typically have mass column densities ~ 0.1 - 0.5 g cm~^ (e.g., 
Arvidsson et al. 2010; Wolf-Chase et al. 2003; Wolf-Chase, Walker, & Lada 1995). All of 
the CARMA cores have mass column densities significantly > 1 g cm~^, the smallest value 
being ~ 6 g cm~^ for core E at Tij = 40 K. 



3.3. CHgOH Masers 



We have discovered three 95-GHz CH3OH masers in the Mol 121 region, coincident 
with MHOs 864, 865, and 892. They are all unresolved by CARMA's 2" resolution, and 
each appear in two contiguous velocity channels. There are very slight (<1") shifts in 
position of the masers between channels, and that, combined with our relatively coarse 
1.538 km s~^ channel width, means that our masers may each be multiple, tightly clustered 
masers, but without better spatial and spectral resolution we will set aside this possibility. 
For each maser, we measured the total flux and estimated the central velocity, and these 
results are presented in Table 5. 

One 44-GHz Class I CH3OH maser was previously discovered in the Mol 121 region (S. 
Kurtz, private communication), at the same location as the brightest 95-GHz maser and 
with a velocity of 2.3 km s"-*^. The channel width of the 44-GHz data was 0.166 km s~^, 
accounting for the discrepancy with our estimated velocity. The 44-GHz entries in Table 5 
were identified from VLA data provided by Kurtz and are discussed in §4.4. 

4. DISCUSSION 

4.1. H2 & [Fe II] Excitation 

Lorenzetti et al. (2002) found that [Fe II] 1.64-//m/H2 2.12-/xm flux ratios for knots 
associated with a jet from IRS 8 in the Vela Molecular Ridge ranged from ~ 0.1 - 1.1 
(uncorrected for extinction). Since the [Fe II] and H2 emitting regions are unsurprisingly 
physically distinct, they computed flux ratios based on two sets of zones, deflned by the 
[Fe II] peak emission and the H2 peak emission. They argue that these high ratios require 
the presence of a fast J-shock component, which produces [Fe II] 1.64-|Um/H2 2.12-/im flux 
ratios in the range 0.5 - 2 for a wide range of parameters (HoUenbach & McKee 1989). Given 
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that the effects of extinction make our observed [Fe II] 1.64-//m/H2 2.12-//m flux ratios lower 
hmits, this suggests the presence of J-shocks in both MHOs 867 & 864; however, the high 
observed H2 2.12-/xm/2.25-/xm ratios across all of the apertures with 2.25-/xm detections 
(6.8 - 14.5) in Table 3 suggest the presence of C-shocks as well. For the high excitation 
temperatures associated with J-shocks, one expects H2 2.12-//m/2.25-//m ~ 3 — 4, and for 
the lower-excitation C-shocks, ~ 6 — 20. Above 20, the excitation temperature is < 1000 
K, and observable emission is not expected (Smith 1994, 1995; Smith et al. 2003). In fact, 
both types of shocks may be found in jets from YSOs. Typically, the molecular emission 
arises from the flanks of a bow shock, and the atomic or ionized component originates close 
to the apex (Smith 2012). 

Neither fluorescence nor shocks can explain a H2 2.12-yum/2.25-/im flux ratio < 1, which 
is the case for the "Anomaly" listed in Table 1 and indicated on Figure 2. We hand-checked 
each image for possible saturation effects near the bright IRAS source. Our combination of 
narrow-band filters and short exposure times (30 s) kept even the brightest sources below 
30k ADU, well within the linear range of the detector (3% linear to 54.4k ADU) and far 
below its full-well capacity of more than 100k ADU. Therefore, we conclude that this effect 
is real. To the best of our knowledge, this is the first time an observed H2 2.12-//m/2.25-//m 
flux ratio < 1 has been reported. Although high-resolution NIR spectroscopy is essential to 
determining the nature of this emission, we note that Dove et al. (1987) found that certain 
circumstances could produce population inversions in para-hydrogen transitions. They 
suggested future work to include the ortho-hydrogen transitions and more completely model 
the processes occurring in shock-heated clouds. Another intriguing possibility is that the 
anomalous emission comes from the immediate front of a fast dissociative shock (T> 10,000 
K), where statistical equilibrium is not reached (M. Smith, private communication). If 
statistical equilibrium is not achieved, the 2.12-yLim/2.25-//m ratio should be the ratio of the 
transition probabilities. The transition probabilities for the H2 2.12-/im 1-0 S(l) and H2 



- 14 - 



2.25-//m 2-1 S(l) transitions are 3.47x10''' s~^ and 4.98x10''' s~\ respectively (Turner et al. 
1977). In this case, the expected ratio would be about 0.7. Although this is still somewhat 
larger than we calculate in Table 1, the ratio is affected by choice of aperture size. The 
Anomaly is bright in both line and continuum NIR images (see Fig. 1), and it is coincident 
with a deeply embedded infrared source (DES) identified from polarimetric patterns in the 
scattered NIR emission (Yao et al. 2000). The DES is also coincident with core C, the 
brightest 3-mm core in the Mol 121 region (see §4.2.3). 

4.2. Dense Cores &: Outflows 

In this section, we discuss each of the five cores we detected with CARMA, previous 
related observations, and the possible association of cores with outflows, individually. There 
is direct evidence of massive star formation occurring in three of the five cores (A, C, & D). 
The remaining two cores (B & E) may be starless. In the absence of significant external 
heating, these cores are hkely to have temperatures < 20 K, in which case the derived 
masses and mass column densities suggest that they may be sites of future massive star 
formation. 

4.2.1. Core A 

Core A is associated with the IRAS source and several signposts of massive star 
formation, including a compact H II region observed at 3.6 cm and 6 cm (Jenness et al. 
1995; Molinari et al. 1998), a PDR detected in our H2 2.12-/xm, 2.25-/im, and [Fe II] 
1.64-/xm observations (Figs. 1 & 2) and evident in the WISE image (Fig. 4), and H2O 
maser emission (Palla et al. 1991; Brand et al. 1994; Anglada et al. 1997). Furthermore, 
core A lies at or near the continuum peak identified through previous sub-millimeter and 
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millimeter observations (FIR 2 - Jenness et al. 1995; McCutcheon et al. 1995; Curran & 
Chrysostomou 2007; Motte et al. 2007). Core A may be associated with a ~E-W outflow 
defined by several MHOs. In particular, the bow-shaped MHO 867 opens toward this 
source, as does MHO 938 further to the west (see Figures 2 & 5). MHO 867 also exhibits 
strong [Fe II] emission. There are several MHO groupings toward the east, which do not 
align in any obvious way (MHOs 939—942). 

4.2.2. Core B 

Core B has not been detected previously with millimeter or sub-millimeter observations, 
nor is it associated with any obvious signposts of star formation such as masers, radio 
continuum, or obvious outflow activity. It is located at the northern end of the bright 
PDR, coincident with an apparent dip in the NIR emission, which suggests high extinction 
towards this core. It may be prestellar in nature. If this is true, its temperature may be 
lower than we assumed, and thus its mass may be larger. 

4.2.3. Core C 

Core C has not been detected previously with millimeter or sub-millimeter continuum 
observations, even though it has the largest integrated flux of the cores we detected with 
CARMA. It does, however, coincide with the location of the DES reported by Yao et al. 
(2000), as well the peak of the anomalous H2 emission. Although Motte et al. (2007) 
mapped the entire Cygnus X molecular cloud complex in 1.2- mm continuum at IRAM, they 
detected only one dense core with their 11" beam within the region we identify as a PDR 
associated with IRAS 20188-1-3928. Our cores A, B, & C all lie within this region. Little 
et al. (1988) flrst detected a bipolar molecular outflow in the vicinity of IRAS 20188-1-3928 
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from CO and HCO+ observations. The molecular outflow is aligned roughly N-S, with the 
redshifted gas to the north. Using a distance of 4 kpc, they noted the outflow properties 
were similar to other massive outflows such as Cep A and NGC 2071. Zhang et al. (2005) 
identified an outflow with the same orientation from CO 2— >1 observations at a similar 
resolution to the observations of Little et al. (1988). They derived significantly smaller 
outflow mass and energetic parameters; however, they used a distance of only 0.3 kpc, which 
is much smaller than the most recent estimate of the distance (Mottram et al. 2011). Zhang 
et al. (2005) noted that the outflow appears to be centered ~6" to the north of the IRAS 
source, coincident with the DES identifled by Yao et al. (2000); however, given their ~30" 
beam and spatial sampling, it is impossible to identify the driving source(s). The image 
of Yao et al. (2000) shows NIR reflection nebulosity oriented generally N-S in the vicinity of 
the IRAS source; however, the polarization vectors exhibit a different pattern in a NW-SE 
direction about the DES (coincident with the bar-shaped extension of the Anomaly evident 
in Figs. 2 & 5). Yao et al. (2000) suggest this may indicate the direction of a newly-burst 
outflow from the DES. The association of core C with the DES identifled by Yao et al. 
(2000), strong millimeter flux, and corresponding lack of cm-continuum emission, make this 
object an excellent candidate MYSO in a pre-UC H II phase of evolution. If this is the case, 
anomalous H2 emission might be a signpost for the early outflow phase of a MYSO, but 
clearly high-resolution NIR spectroscopy is necessary to search for any such evolutionary 
trends. 

4.24. Core D 

Core D is associated with a H2O maser and faint continuum emission at 3.6 cm and 
6 cm, indicative of the presence of a UC H II region (Jenness et al. 1995; Mohnari et 
al. 1998). Using the UKT14 receiver at the JCMT, Jenness et al. (1995) detected a 
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submillimeter source (FIR 1) at a position between our cores D & E from 450 //m & 800 //m 
~ 18"-resolution observations. MHOs 864 & 865, the brightest MHOs in the region (with 
the exception of MHO 867, which may be contaminated with fluorescent emission from 
the PDR), exhibit bipolar morphology in a NE-SW direction about core D. Both MHOs 
are associated with CH3OH maser emission at 95 GHz & 44 GHz. We have also detected 
[Fe II] emission at two positions within MHO 864. Several MHOs lie near or along a slightly 
different position angle from the direction defined by MHOs 864 & 865: to the NE, MHOs 
863, 891 & 892, and to the SW, MHO 937. MHO 892 is also coincident with CH3OH maser 
emission at 95 GHz & 44 GHz. Given the presence of both radio continuum emission and 
one or more well-coUimated H2 outflow, this object may be close to transitioning to a stage 
where an expanding UC H II region destroys the MYSO(s)' accretion disk. The bright 
MHO 866 lies directly to the NW of core D, but its association is unclear. 

4.2.5. Core E 

The position of core E coincides with the position of a dense core detected by Motte et 
al. (2007) in 11 "-resolution 1.2-mm continuum observations using the IRAM 30-m telescope. 
For an assumed distance of 1.7 kpc and temperature of 20 K, Motte et al. (2007) derive a 
core mass of 40 M© and density, n(H2), of 3.5x10^ cm~^. For the same temperature, we 
derive a core mass of 12.8 Mq. The difference can be explained if the CARMA observations 
resolve out some of the large-scale structure revealed in sub-millimeter JCMT maps 
(McCutcheon et al. 1995; Jenness et al. 1995; Curran & Chrysostomou 2007). MHOs 
893—897 lie in the vicinity of core E, but it is not possible to associate any of these with a 
driving source presently. 



4.3. H2 Luminosity 



WSS12 noted that the MYSO, Mol 160, exhibits a very high hnjl^boi ratio, faUing 
approximately a factor of two above the relationship that was determined for low-mass 
YSOs (Caratti o Garatti et al. 2006) and extended to MYSOs based on results for IRAS 
20126+4104 (Caratti o Garatti et al. 2008). Although it is not generally possible to 
link individual outflows with driving sources in massive star-forming regions, and caution 
must be exercised in attempting to draw conclusions about the nature of MYSO outflows 
from the Lj/j/Lfeoi relationship alone (WSS12), it is instructive to compare the total H2 
luminosity from all detected outflows to the total bolometric luminosity of the Mol 121 
region. If the H2 flows are produced by low- or intermediate-mass YSOs, one would expect 
the H2 emission to be under-luminous with respect to the bolometric luminosity, which 
comes primarily from the IRAS source. 

There are many uncertainties in estimating the total H2 luminosity, including the 
distance to Mol 121 /IRAS 20188-1-3928, the extinction toward individual H2 knots, and the 
ratio of total H2 to H2.12 luminosity. Nevertheless it is possible to place reasonable limits 
on the assumptions. Without spectra of the MHOs, we can only estimate the effects of 
extinction. Although the extinction is very large in the dense cores, most of the MHOs 
are located in the more extended cloud. It seems reasonable to assume the column density 
derived for the parsec-scale NH3 clump in this region by Anglada et al. (1997), N(H2) = 
5.5x10^^ cm~^. Using the relation Ay/N(H2) ~ 5.3 xlO~^^ cm^ mag yields an estimate 
of the visual extinction Ay fa 29 mag. Further assuming an optical parameter for dense 
clouds of Ry = 5, yields A2.12 ~ 3.6 mag (Allen 2000). 

The luminosity of the H2 2.12 //m emission is given by: 

L2.12 = 47rd'F2.i2 X 10°-^^2-i2 (3) 
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F2.12 = 4.59 X 10"^^ W m~^, summed over all the MHOs listed in Table 1 (including the 
PDR and Anomaly increases F2.12 by a factor of ~ 2.5.) Depending upon the temperature 
of the shocked H2, the 2.12-/xm luminosity is typically ~ 5-10% of the total ro-vibrational 
H2 emission, Lh^ (e.g., Caratti o Garatti 2006, 2008). This yields Lh^ ~ 11-7 - 23.3 for 
the combined MHO luminosity. The bolometric luminosity of IRAS 20188+3928, calculated 
from the far- infrared SED and assuming d = 1.7 kpc, is 8.2 xlO^ Lq (Mottram et al. 
2011). These results place Mol 121 a factor of 1.5 to three above the Li^j/Lfcoi relationship 
determined by Caratti o Garatti et al. (2006, 2008). This is very similar to the result 
obtained by WSS12 for Mol 160, and strongly suggests that the H2 flows are driven by 
massive YSOs. 

4.4. CH3OH Masers 

Class I CH3OH masers are well-correlated with molecular outflows in massive 
star-forming regions (Cyganowski et al. 2009 and references therein; Fontani et al. 2010) 
and with the presence of compact 3-mm continuum emission (Schnee & Carpenter 2009), 
suggesting that these masers are signposts of an early stage in the evolution of a MYSO 
before an expanding UC H II region has destroyed the accretion disk. 

The strongest CH3OH maser detection (Ml) in Mol 121 is associated with MHO 864 
(Figure 5) . It is by far the brightest at 95 GHz and was the only previously detected maser 
at 44 GHz (S. Kurtz, private communication). MHO 865 and the associated maser (M3) 
are positioned exactly opposite core D, while MHO 892 and its associated maser (M2) lie 
along a slightly different position angle, though they also appear to be associated with 
core D. The masers all have low velocities, close to the Vlsr of the ambient cloud, which 
is reported to be 1.5 km s"^ by Zhang et al. (2005) and listed as 2.1 km s~^ in the RMS 
catalog (Mottram et al. 2011). This is consistent with the suggestion that CH3OH masers 
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arise in systemic gas in outflow cavity waUs (DeBuizer et al. 2009). 

WSS12 suggested that the 95-GHz/44-GHz ratio might prove to be a useful diagnostic 
for the shortest-hved (e.g., very early) phases of massive star formation. In order to 
estimate the flux ratios of our masers, we binned the 44-GHz data (provided by S. Kurtz) 
to approximate our wider velocity channels. This summed data lowered the noise and 
revealed 44-GHz maser emission at the location of all three of our masers, not just the one 
known previously. We were then able to measure and compare the fluxes at 44 GHz and 95 
GHz. These fluxes are all presented in Table 5. All three masers arc brighter at 95 GHz 
than at 44 GHz, in contrast to the statistical result reported by Fontani et al. (2010), who 
concluded that the 95-GHz line is intrinsically fainter based on their detection rates. Ml 
and M3, which are associated with the same outflow angle, have remarkably similar flux 
ratios (3.73 and 3.69, respectively). M2, which is associated with MHO 892, has a 95-GHz 
to 44-GHz flux ratio of 6.63, dramatically higher than the other two. MHO 892 is much 
more compact than MHOs 864 & 865, though we do not know if this is related in any way. 

4.5. Different Stages of Massive Steir Formation in Mol 121 

We can tentatively order the dense cores in Mol 121 by evolutionary stage, based 
on their continuum and spectral-line properties. (1) Core A is associated with the IRAS 
source, a compact PDR, and bright cm-continuum emission, suggesting it is the oldest site 
of massive star formation in Mol 121. (2) Core D, with its associated faint cm-continuum 
emission, luminous H2 outflow, CH3OH masers, and Fe II emission may harbor one or more 
massive objects that has not yet destroyed its accretion disk. (3) Core C, with its DES, 
likely contains a MYSO in a pre-UC H II region phase of evolution. (4) B & E appear to be 
massive starless or prestellar cores. We note that cores B & C also lie near the periphery 
of the PDR, which suggests that Mol 121 would be an excellent candidate for follow-up 
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studies to explore the possibility of star formation triggered by the expanding H II region 
associated with IRAS 20188+3928. 

5. SUMMARY AND CONCLUSIONS 

1. Using NICFPS on the ARC 3.5-m telescope at APO, we have discovered 12 new MHOs 
in the vicinity of Mol 121/ IRAS 20188+3928, bringing the total known MHOs in this 
region to 17. We have been able to separate one of the five previously discovered 
MHOs from a FDR and a region marked by an anomalous H2 2.12-//m/2.25-//m flux 

ratio < 1. To the best of our knowledge, this is the first time an observation of a H2 
2.12-/im/2.25-/xm flux ratio < 1 has been reported. 

2. We have detected [Fe II] 1.64-/im emission toward MHOs 864 & 867. [Fe II] 
1.64-/xm/H2 2.12-/xm and H2 2.12-/xm/2.25-/xm flux ratios measured along and across 
these MHOs indicate the presence of both J-shock and C-shock components. 

3. CARMA 3-mm continuum and line observations have revealed the presence of flve 
compact cores (A, B, C, D, & E) and three CH3OH masers (Ml, M2, & M3) in the 
Mol 121//i?^5' 20188+3928 region. The physical properties of the cores indicate that 
cores A, C, & D are presently forming massive stars, while cores B & E may be the 
sites of future massive star formation. Masers Ml, M2, & M3 he toward MHOs 864, 
892, & 865; all of which are associated with outflow activity from one or more sources 
embedded in core D. The 95-GHz masers all have 44-GHz counterparts, but they are 
all brighter at 95 GHz. The two masers (Ml & M3) associated with bipolar MHOs 
864 & 865 have strikingly similar ratios of ~ 3.7, while M2 has a higher ratio of ~ 
6.6. The significance of this is presently unknown, but suggests that future studies 
of maser line ratios are warranted, and may hold important implications for probing 
physical conditions in massive star-forming regions. 
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4. The positions and orientations of the MHOs relative to the cores suggest the presence 
of multiple outflows originating from at least two of the cores. A comparison of the 
total H2 luminosity of the MHOs to the total bolometric luminosity of this region 
strongly suggests the outflows are driven by massive YSOs. 

5. We tentatively age-order the cores from most to least evolved based on the different 
signposts of star- formation activity. (1) Core A is associated with the IRAS source, a 
compact PDR, and bright cm-continuum emission, suggesting it is the oldest site of 
massive star formation in Mol 121. (2) Core D, with its associated faint cm-continuum 
emission, luminous well-coUimated H2 outflow(s), CH3OH masers, and Fe II emission 
may harbor one or more massive objects that has not yet destroyed its accretion disk. 
(3) Core C is not associated with cm-continuum emission, but it coincides with a 
DES, the anomalous H2 2.12-/xm/2.25-/xm emission, and possible outflow. It likely 
harbors a MYSO in a pre-UC H II region phase of evolution. (4) While cores B & E 
are massive, they also appear to be starless or prestellar. 
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Table 2. Mol 121 [Fe II] Emission Line Fluxes 



Associated MHO Area a(J2000)'' 5(J2000)* Fi,e4^rn ^ 

10-10 Sterad (h m s) (° ' ") lO"!* W m'^ 



864 1.833 20 20 39.8 +39 38 22 7.04±1.17 

0.668 20 20 39.5 +39 38 14 2.14+0.36 

867 3.564 20 20 38.2 +39 37 51 19.89+4.07 

^Coordinates correspond to emission peaks. 

''Flux errors combine three categories of uncertainty (calibration, subtraction, 
and aperture) and correspond to +1 cr (see §2.1). 
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Table 3. MHO 864 & MHO 867 Emission Line Flux Ratios 



MHO # Aperture # a(J2000)=' (5(J2000)'' F1.64/F2.12 F2.12/F2.25 
(h m s) (° ' ") 



1 


20 


20 


39.8 


+39 


38 


23 


0.9+0.2 


>5.9 


2 


20 


20 


39.8 


+39 


38 


21 


0.8+0.1 


7.0+1.8 


3 


20 


20 


39.7 


+39 


38 


19 


0.1+0.03 


6.8+1.3 


4 


20 


20 


39.7 


+39 


38 


17 


<0.03 


13.8+2.8 


5 


20 


20 


39.6 


+39 


38 


15 


<0.03 


14.5+4.5 


6 


20 


20 


39.6 


+39 


38 


14 


0.4+0.1 


12.8+4.3 


7 


20 


20 


39.8 


+39 


38 


17 


<0.04 


>7.9 


8 


20 


20 


40.0 


+39 


38 


17 


<0.3 


>1.4 


9 


20 


20 


39.6 


+39 


38 


19 


<0.03 


>7.1 


10 


20 


20 


39.4 


+39 


38 


19 


... b 


>0.7 


11 


20 


20 


38.2 


+39 


37 


55 


1.2+0.3 


>4.0 


12 


20 


20 


38.1 


+39 


37 


53 


0.5+0.1 


>7.2 


13 


20 


20 


38.1 


+39 


37 


52 


0.6+0.1 


10.4+3.6 


14 


20 


20 


38.2 


+39 


37 


50 


0.5+0.1 


>6.5 


15 


20 


20 


38.3 


+39 


37 


53 


1.2+0.2 


>4.5 


16 


20 


20 


38.2 


+39 


37 


51 


1.8+0.3 


>5.6 


17 


20 


20 


38.4 


+39 


37 


52 


1.9+0.5 


>2.3 



''Coordinates of aperture centers. 
''ct < 3 
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Fig. 1. — Inverted greyscale images of H2 2.12-//m (top left), H2 2.25-//m (top middle), and 
[Fe II] 1.64-/xm emission (top right); and continuum-subtracted H2 2.12-/im (bottom left), 
H2 2.25-/xm (bottom middle), and [Fe II] 1.64-/xm (bottom right) emission in the Mol 121 
region. The cyan X marks the position of the IRAS source. MHOs reported by Varricatt et 
al. (2010) are labelled and enclosed in red and MHOs reported in this paper are labelled and 
enclosed in green on the continuum-subtracted H2 2.12-|Um image. [Fe II] knots are enclosed 
with blue dashed lines on the continuum-subtracted [Fe II] 1.64-/xm image. 

Fig. 2. — Inverted greyscale image of continuum-subtracted H2 2.12-/xm emission. MHOs 
reported by Varricatt et al. (2010) are labelled and enclosed in red, MHOs reported in 
this paper are labelled and enclosed in green, [Fe II] knots are enclosed with blue dashed 
lines, and peaks of 3-mm cores identified with CARMA are indicated with magenta crosses 
and labelled as discussed in the text. The PDR and anomalous H2 emission (Anomaly) are 
labelled. 

Fig. 3. — Inverted greyscale image of continuum-subtracted H2 2.12-//m emission in the 
vicinity of MHO 864 (left) and MHO 867 (right). MHOs and [Fe II] knots are enclosed as in 
Figure 2. Aperture regions corresponding to entries in Table 3 are indicated and labelled. 

Fig. 4. — Three-color WISE image of the Mol 121 region. Colors are 3.4 fim (blue), 4.6 fim 
(green), and 12 /im (red). MHOs are indicated. The black X denotes the position of the 
IRAS source. 

Fig. 5. — Inverted greyscale of continuum-subtracted H2 2.12-/xm sub-image. CARMA 3- 
mm continuum cores (labelled green contours), [Fe II] emission peaks (yellow circles), and 
95-GHz CH3OH masers (labelled cyan Xs) are indicated. Red crosses show the positions of 
VLA 3.6-cm sources reported by Jenness et al. (1995). 



